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Summary 

Numerical analysis and experimental results 
are presented to define a method for quantitatively 
measuring the temperature distribution of a spheri- 
cal diffusion flame using Rainbow Schlieren 
Deflectometry in microgravity. First, a numerical 
analysis is completed to show the method can 
suitably determine temperature in the presence of 
spatially varying species composition. Also, a 
numerical forward-backward inversion calculation 
is presented to illustrate the types of calculations 
and deflections to be encountered. Lastly, a nor- 
mal gravity demonstration of temperature meas- 
urement in an axisymmetric laminar, diffusion 
flame using Rainbow Schlieren deflectometry is 
presented. The method employed in this paper 
illustrates the necessary steps for the preliminary 
design of a Schlieren system. The largest deflec- 
tions for the normal gravity flame considered in 
this paper are 7.4x10 radians which can be 
accurately measured with 2 m focal length colli- 
mating and decollimating optics. The experimental 
uncertainty of deflection is less than 5xl0 -5 
radians. 


Introduction 

Optical diagnostic techniques can be classified 
as either local or path integrated. Local measure- 
ment techniques for microgravity combustion 
research include (1) Laser Doppler Velocimetry 
(LDV) and Particle Image Velocimetry (PIV) to 
measure the velocity field, and (2) Rayleigh scat- 
tering and Laser-Induced Fluorescence (LIF) or 
Planar Laser-Induced Fluorescence (PLIF) to 
measure temperature and/or species concentrations 


(ref. 1). Application of these techniques presents 
unique challenges and difficulties in the micro- 
gravity environment where mass, volume, and 
electrical power available are highly limited and 
severely constrained. In addition, remote operation 
is also often required further complicating the 
implementation. The LDV method enables de- 
tailed velocity information at a point but in order 
to determine the entire velocity field Particle 
Image Velocimetry (PIV) (ref. 2) is often applied. 
However, PIV methods are typically not well 
suited to characterize extreme turbulence because 
of low spatial and temporal resolution and cannot 
account for large out of plane motions. Seeding of 
the flow is a challenge in microgravity as well. 
The scattering and fluorescence methods require 
bulky, power consuming lasers often unsuitable 
for use in microgravity, and are sensitive to mis- 
alignment. These methods will become more 
applicable to microgravity as technology advances 
in solid-state lasers and advanced imaging 
systems. 

The path-integrated techniques used for 
microgravity overcome some of the drawbacks of 
the local measurements techniques. The measure- 
ments must be inverted tomographically to obtain 
local distributions of scalar or flow properties of 
interest. The inversion process is relatively simple 
for laminar flows if symmetry can be assumed. 
But for turbulent flows the deconvolution is tedi- 
ous. Examples of path-integrated methods applied 
in microgravity include emission/ absorption 
imaging pyrometry (ref. 1), for measurement of 
temperature, CH emission tomography (refs. 3 
and 4), infrared diode laser wavelength modula- 
tion spectroscopy (WMS) for measurement of 
oxygen and water vapor concentrations (refs. 5 
to 8), shearing plate interferometry (ref. 9), Fan- 
Beam Emission Tomography (FBET) (ref. 10), 
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and quantitative rainbow Schlieren deflectometry 
(ref. 11). The objective of this paper is to discuss 
the application of quantitative rainbow Schlieren 
deflectometry for the case of a laminar spherical 
flame. 


Experiment Description 

The Schlieren method is one of the simplest and 
oldest methods for visualizing refractive-index 
inhomogeneities in transparent media with its 
origins dating back to the 17th century. The 
Schlieren instrument is sensitive to transverse 
refractive-index gradients in a media. The present 
application is a laminar flame. The gradients cause 
angular deflections of the incoming light rays 
which are transformed at the filter plane of the 
Schlieren system by a spatial filter. Many different 
types of spatial filters have been developed to 
relate the angular deflections to different quanti- 
ties. The various types of spatial filters, Schlieren 
optical configurations and applications are exten- 
sively reported in the literature (refs. 12 to 15). 

A simple Lens-Type Schlieren System is 
shown in figure 1. An input lens collimates the 
light from a light source through a transparent 
media with transverse gradients of index of refrac- 
tion, referred to as the test section. The spherical 
flame is located in the test section. As the light 
exits the test section, it is angularly deflected 
according to a line integral eq. (1). The optical 
inhomo genities refract or bend light rays in proportion 
to their gradients of refractive index. 

Q exll = tan Q exit = j>V/?^ (1) 

L 

A decollimating lens refocuses the light onto a 
filter plane. A continuous rainbow filter shown in 
figure 2, first introduced by Howes (ref. 16) alle- 
viates undesirable optical discontinuities, and 
avoids the nonlinear sensitivity associated with 
conventional grey-scale methods. An image of the 
test section and the rainbow filter is then focused 
by a camera lens onto the color sensitive detector, 
allowing accurate quantification for the color 
attributes of the resulting image, and hence, the 
associated ray deflections (ref. 1 1). The “color” of 
the resulting Schlieren image is quantified using 
an RGB relationship as measured using the image 



Figure 1. — Schematic drawing of a simple lens-type 
Schlieren system. 



Figure 2. — Typical circular Schlieren rainbow color 
filter used in microgravity experiments for spherical 
flames; filter diameter is 8 mm and “clear” center 
diameter is 0.5 mm. 


Diffusion flame 


t 

Fuel (oxidizer) 

Figure 3. — Schematic representation of the spherical 
burner generated gaseous diffusion flame experiment. 

hue (ref. 11). The spherical burner arrangement is 
schematically represented in figure 3. Typically, a 
porous sphere is attached to a “fuel” supply tube 
and fuel is uniformly dispersed in the vicinity of 
the burner. A quiescent ambient “air” or other 
oxidizer diffuses to the reaction zone of the diffu- 
sion flame. In microgravity the buoyant effects are 
significantly eliminated thereby allowing a spheri- 
cally symmetric, one -dimensional formulation. 
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Numerical Analysis to Demonstrate 
Temperature Sensitivity 

In order to investigate the optical properties of 
the flame, both the gas composition and tempera- 
ture must be known (ref. 17) since both affect the 
index of refraction. A time-dependant, one- 
dimensional numerical model based on Chemkin 
was used to determine flame temperature and 
species concentrations (ref. 18). Figure 4 shows a 
computed radial temperature profile with this 
model after 5.4 sec from ignition with a fuel 
composition of 20% H2/25% CH4/55% N 2 and a 
fuel flow rate of 24 mg/s. This flame is burning in 
pure air with 21 percent oxygen and 79 percent 
nitrogen. The hottest part of the flame, 2000 K, is 
located at 5 cm from the center of the spherical 
flame. The ambient air is constant at 300 K from 
13 to 20 cm indicating that no temperature gradi- 
ents are present at this position in the domain. The 
mole fractions of the major species in the flame 
are nitrogen, oxygen, methane, and hydrogen, are 
shown plotted in figure 5 over a 20 cm radial 
domain after 5.4 sec from ignition. The fuel, 
methane and hydrogen, and oxidizer, oxygen, are 
completely consumed at a radial position of 5 cm. 
The fuel rich side of the flame occurs over the 
radial domain of 0 < r < 5 cm. The oxidizer rich 
side of the flame occurs when r > 5 cm. At radial 
locations greater than 1 3 cm no gradients in flame 
“composition” are observed. The radial profiles of 
other important major species in the flame are 
shown in figure 6. Water vapor, carbon dioxide, 
and carbon monoxide are produced in the reaction 
zone of the diffusion zone of the flame. Carbon 
monoxide is generated on the fuel rich side of the 
flame but is consumed in the high temperature 
reaction zone at radial locations approaching 5 cm. 
The peak values of water vapor and carbon diox- 
ide gas are produced in the highest temperature 
regions of the flame near 5 cm and are observed to 
diffuse toward both the fuel and oxidizer sides of 
the diffusion flame. 

In the determination of the index of refraction 
in the flame an analysis in which all of the com- 
puted major and minor species included in table I 
have been considered in the calculation of the 
index of the refraction profile. This analysis is 
used to answer the following two important 
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Figure 4. — Computed temperature profile with time 
dependent Chemkin model for a spherical diffusion 
flame burning in air at 1 .2 atm and 300 K; fuel type: 
20% H 2 / 25% CH 4 / 55% N 2 with mass flow rate of 
24 mg/sec, time from ignition is 5.4 sec. 
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Figure 5. — Fuel and oxidizer specie mole fractions for 
an adiabatic calculation with fuel 20% H 2 / 25% CH 4 / 
55% N 2 in air, 1 .2 atm, fuel mass flow rate 24 mg/sec, 
at 5.4 sec after ignition. 


18x10-2 
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Figure 6 . — H 2 O, CO 2 and CO species mole fractions for 
an adiabatic calculation with fuel 20% H 2 / 25% CH 4 / 
55% N 2 in air at 1 .2 atm, fuel mass flow rate 
24 mg/sec, at 5.4 sec after ignition. 
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TABLE I.— INDEX OF REFRACTION AT 273.15 K, 

1 ATM FOR VARIOUS GASES OF INTEREST (REF. 19). 


Values are relative to a vacuum and for a temperature of 
273.15 K, wavelength 

A, = 590 nm and 760 mm pressure 

Specie 

Formula 

Index of refraction, 
n 

Air 

Air 

1.000293 

Carbon Dioxide 

co 2 

1.000451 

Carbon Monoxide 

CO 

1.000338 

Helium 

HE 

1.00036 

Hydrogen 

h 2 

1.000132 

Hydrogen Atom 

H 

1.000407 

Methane 

ch 4 

1.000444 

Methanol 

CH 3 OH 

1.000586 

Nitrogen 

N 2 

1.000297 

Oxygen 

0 2 

1.000272 

Oxygen Atom 

0 

1.000649 

Water 

h 2 o 

1.000254 


questions concerning the sensitivity of Schlieren 
deflectometry to temperature and gas composition 
(ref. 17). 

(1) How do chemical species affect the index 
of refraction? 

N 

<n(T)-l) mix = £ Xi * (n i (T)-\) i (2) 

i=l 

(2) How does temperature affect the index of 
refraction? 

("- l )ho,=( n o-')coid*y ( 3 ) 

Where n is the index of refraction, X, are the mole 
fraction of specie i, N is the total number of gas 
species, n a is the reference index of refraction 
(ref. 19), listed in table I, T 0 is the temperature used 
for the referenced index of refraction, and T is the 
local gas temperature. Other useful data for the 
refractivity of combustion gases can be found in 
ref. 25. Knowing the relative contribution of 
chemical composition and the temperature on the 
refractive-index is necessary for quantitative flame 
measurements. 



Radius (cm) 


Figure 7. — Index of refraction n(r) profile of the flame 
with complete species and air as calculated with the 
Chemkin computed species and temperature profiles, 
20% H 2 / 25% CH 4 / 55% N 2 in air at 1 .2 atm, fuel mass 
flow rate 24 mg/sec, at time from ignition of 5.4 sec. 

Using the data in figures 4 to 6 and eqs. (2) 
and (3), the refractive index radial profile can be 
computed. 

Figure 7 shows a comparison between the 
refractive index for air and all of species shown in 
table 1, with the same computed temperature 
profile. The analysis was conducted on several 
profiles with similar results represented in 
figure 7. Figure 7 shows that the index of refrac- 
tion for the flame is nearly equivalent to high 
temperature air and therefore the data reduction of 
refractive-index to temperature can be made to 
suitable accuracy by neglecting variations in 
composition. This illustrates the significant effect 
of temperature rather than gas composition on the 
variation of index of refraction through the flame. 
The reasoning for this can be understood by noting 
the flame is composed of primarily nitrogen, 
which is also a major component in air. Thus, the 
gas temperature has the dominant effect on refrac- 
tive index rather than the gas composition. 

Forward-Backward Inversion 
Calculation on Synthetic Data 

The tomographic reconstruction of an axi- 
symmetric refractive-index distribution from 
angular deflections (e.g., Schlieren data) for a 
spherical diffusion flame is demonstrated with a 
synthetic data set. The purpose of this exercise is 
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to show how temperature can be extracted from 
deflection data as measured from Rainbow 
Schlieren deflections. In addition, this exercise 
assists in the design of the Schlieren system by 
determining both angular and linear deflections. 
This is conducted for the microgravity spherical 
diffusion flame described in the previous section 
and shown schematically in figure 3. 

The angular deflections generated by the 
spherical diffusion flame can be computed from 
the radial refractive-index distribution shown in 
figure 7 using the forward expression (ref. 20). 


e(y)=2 yj 


d(n — l) dr 


dr 


2 2 
r ~y 


12 


(4) 


6X10- 4 



Radial position (cm) 

Figure 8. — Angular defection of light rays through the 
flame shown in figures 4 to 7 as computed using 
equation (4). 


The corresponding transverse deflections 
occurring in the plane of the rainbow filter can be 
computed for a known focal length decollimating 
lens (ref. 11). 


d f (y) = F L E(y) (5) 

These calculations have been completed in a 
MATHCAD 2000 program worksheet and the 
improper integral in eq. (4) has been evaluated 
with an infinite limit approximation. 

Figure 8 shows angular deflection of light rays 
through the microgravity spherical flame. Over the 
radial position of 0 < /- < 3 cm the deflections are 
negative and the flame acts like a “focusing” lens. 
For radial stations of 3 cm < r < 12 cm the flame 
disperses the light rays similar to a defocusing 
lens. The maximum angular deflection is ap- 
proximately 5.4X10" 4 radians and appears at 9 cm 
from the center of the flame. If the ambient condi- 
tions in this flame were constant at a radius of 
20 cm the angular deflection should be zero. 
Examination of figure 8 shows that the angular 
deflections are in fact nonzero and negative at 
20 cm. Physically, however, these deflections 
should be zero since the flame has constant tem- 
perature and gas composition, i.e., pure air, at a 
radius of 20 cm as observed in figures 4 to 6. The 
reason for this discrepancy is due to curve fitting 
of the index of refraction profile, figure 7, which 
was required in order to make use of the 
MATHCAD 2000 program. As a result, a 


numerical error appears in figure 8 for a radius 
greater than 15 cm. This curve fitting error, 
although not physically correct, does not severely 
impact the final results appearing in figures 10 and 
1 1 . This is actually an interesting result since this 
demonstrates that experimental uncertainties in 
measured angular deflections in the ambient gas 
do not result in large errors of temperature. Also, 
this shows that the field of view of the Schlieren 
image need not have constant index of refraction 
in the ambient condition in order to obtain a 
quantitative result. 

If a decollimating lens of focal length 3 m is 
assumed, eq. (5) can be applied to the calculated 
angular deflections shown in figure 8 to obtain the 
linear deflections at the filter plane. Figure 9 
shows these linear deflections at the filter plane 
enabling sizing of the color filter. The maximum 
linear deflection is approximately 1.7 mm so a 
rainbow filter of at least 4 mm in diameter would be 
required for this application. Since this flame is 
spherically symmetric an Abel-Transform inversion 
can be performed to reconstruct the refractive-index 
field from the deflection information shown in 
figure 8. The Abel Transform for Refractivity, 
8(r) = n(r) - 1 (ref. 20), is given as: 



where 8(r) = n(r)-l 


(6) 
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Figure 9. — Linear deflection at the filter plane for the 
flame shown in figures 4 to 7 as computed using 
equation (5) for a 3 m focal length decollimating optic. 



Figure 10. — Comparison of numerically generated index 
of refraction field based on the Chemkin model with 
the Abel transformed deflections from figure 7 (eq. (6)) 
for the flame shown in figures 4 to 7. 


Figure 9 shows the refractive-index field computed 
using an Abel Transform Inversion calculation 
using eq. (6) (red line) and the refractive -index 
calculated from flame properties and shown in 
figure 7 (blue line). The observed agreement is 
very good, demonstrating the numerical calcula- 
tions performed and the deflections computed in 
figures 8 and 9 are correct. The gas temperature 
can be rapidly computed from figure 10 using 
eq. (3) assuming the gas composition is air. As can 
be seen in figure 11, the Chemkin generated 
temperature profile in figure 4 (blue line) and the 
temperature profile derived from the backward- 
inversion (i.e., Abel Transform) calculation of the 
synthetic data (red line) produce nearly the same 
result, demonstrating the methodology. 

This exercise also serves to provide a prelimi- 
nary design for a Rainbow Schlieren system. The 
maximum angular deflection for a spherical diffusion 
flame is approximately 5.4X10" 4 radians, such that a 
decollimating lens having a focal length of 3 m 
provides a corresponding maximum lateral deflec- 
tion in the plane of the rainbow filter of 1.7 mm. 
Filters of similar dimensions have been previously 
fabricated. As further verification of this calculation, 
previous quantitative temperatures have been meas- 
ured in the 2.2-second drop tower at NASA Glenn 
Research Center using Rainbow Schlieren deflecto- 
metry in flames. In these experiments (refs. 21 and 
22) a low Reynolds number, pure hydrogen fuel gas 
jet was allowed to bum in microgravity. The maxi- 
mum angular deflections measured were approxi- 
mately 5X10” 4 radians, in reasonable agreement with 
the calculation shown in figure 8. 



Radius (cm) 


Figure 1 1 . — Comparison of the flame temperature 
profile with Chemkin model (fig. 4) and Abel 
transformed angular deflection (fig. 8). 


Demonstration Experiment in 
Normal Gravity 

A normal gravity demonstration experiment 
was completed in a steady air coflow flame under 
atmospheric, normal gravity conditions. The fuel 
used was 20% 1^25% CtU/55% N 2 at 9.3 ml/s in 
a central fuel tube with a surrounding coflowing 
air at 583 ml/s. The coflow flame is open to the 
atmosphere and is burning in quiescent, ambient 
air. The burner specifications have been designed 
to match the burner used by Santoro et al. 
(ref. 23). Figure 12 shows the blue flame with a 
visible flame height of approximately 35 mm and 
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e(r): Angular deflection (radians) 



Figure 1 2. — Visible flame image of “1 -g” co-flow flame 
used for Schlieren demonstration experiment for 
20% H 2 / 25% CH 4 / 55% N 2 : fuel flow rate: 9.3 ml/sec 
air flow rate: 583 ml/sec. 



Figure 13. — Color Schlieren image of flame shown in 
figure 1 2 with a symmetric color filter. 



Distance from centerline (mm) 


Figure 14. — Measured angular deflections of light rays 
passing taken directly from figure 13 at 9 mm above 
the burner base. 


flame base of 12 mm. A color Schlieren image 
was collected with a symmetric filter, shown in 
figure 13 with a Z-pattem Schlieren system with 
two off-axis spherical mirrors of focal length 2 m. 
A symmetric, noncircular rainbow filter of 2 mm 
half-width was used to obtain the Schlieren image 
in figure 13. The “red” background has zero 
deflection. The measured angular deflection data 
across the flame at 9 mm above the flame base is 
shown in figure 14. The maximum deflection is 
approximately 7x1 (T 4 radians, or 1.4 mm on the 
filter, at a radial position of 10 mm from the 
central axis of the flame. The Schlieren deflection 
data shows that “perfect” symmetry is not com- 
pletely achieved but is a reasonable assumption. 
The experimental or instrument uncertainty is less 
than 5xl0 -5 radians or 0.1 mm linear deflection, 
which is very accurate. The angular deflection data 
shown in figure 14 was then transformed using a 
simple approximation to the integral shown in 
eq. (6) (ref. 24). 

Figure 15 shows the thermocouple measure- 
ments for comparison with the Schlieren at a 
9 mm height above the burner base. The thermo- 
couple measurements have been corrected for 
radiation losses which are 100 K at most in the 
hottest part of the flame located at a radial position 
of 6 mm. It is noteworthy that the maximum 
deflection at r = 10 mm is not equal to the position 
of maximum flame temperature at r = 6 mm 
similar to the numerical calculations presented 
earlier. This is true since Schlieren deflectometry 
responds to gradients in the refractive index. 
A comparison of the thermocouple data and the 
Schlieren results are presented in figure 16. The 
results in figure 16 show reasonable agreement 
with measured thermocouple data and the 
Schlieren determined temperatures accept near the 
center of flame. Improvements can be realized in 
the comparison by improvements in the algorithm 
to approximate the numerical integration of 
eq. (6). However, the largest errors will always 
occur near the center line since the integral in 
eq. (6) has a mathematical singularity at the radial 
position of zero. In other words, the solution of 
eq. (6) for Refractivity blows up as r approaches 
zero (i.e., as r — > 0, 5 — > <») so results at small radii 
will generally over predict the actual value. 
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T: Temperature (K) ti T: Temperature (K) 
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gure 1 5. — Radiation corrected flame temperatures measured with thermocouple in flame shown 
n figure 1 2 at 9 mm above burner base. 
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Figure 16. — Temperature comparison between Abel transformed Schlieren data and thermocouple 
measurements of the flame shown in figure 12 at 9 mm above the burner. 
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Conclusions 

The following conclusions can be made from 
the results presented. 

1. For the case of a fuel blend 20% H 2 / 25 % 
CH4/55% N 2 , the gas can be treated ana- 
lytically equivalent to heated air when per- 
forming the Abel transformation since the 
flame is composed of mostly nitrogen. For 
other fuels, if the major species concentra- 
tions are known, eq. (2) and (3) enable 
accurate determination of the flame 
temperature. 

2. The methods used in this paper enable a 
designer to determine the primary parame- 
ters for a Schlieren optical system, (a) angu- 
lar deflections (b) focal lengths of the 
collimating and decollimating lenses or mir- 
rors, and (c) size of the rainbow filter, based 
upon flame properties. A numerical model 
for chemical species and temperature 
(ref. 18) as well as and forward/backward 
Abel inversion calculation enables this. 

3. A normal gravity demonstration experiment 
was successfully completed to illustrate the 
methodology of Quantitative Rainbow De- 
flectometry for a flame. A comparison of the 
temperature profile in the flame between a 
radiation corrected thermocouple and the 
Abel inverted Schlieren data are in reason- 
able agreement. 
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